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1 Introduction 

In the first years after tlie discovery of superfluid He-4, tlie property of superfluidity attracted 
attention not only due to its uncommonness but also to the hope of some physicists for its 
connection with the basically new, still undiscovered laws of the microworld [1]. After seven 
decades from the discovery of the superfluidity (P. Kapitsa [5], J- Allen, A. Misener |3]) and a 
century from the time when He-4 was fluidized (H. Kamerlingh-Onnes [1]), it becomes clear 
that superfluid helium justified hopes for the manifestation of a number of unusual properties, 
but no "new physics" has arisen because of helium. 

In what follows, I will characterize briefly the state of the theory of superfluid IIe-4 (He- 
ll) and mark the unsolved questions. In this case, I do not pretend to the completeness of 
analysis or references and do not try to avoid the own sympathies for some questions. The 
other reviews on the theory of superfluidity He-II are, for example, [5]-[7], and a short one is 
presented in [S]. 

2 History 

A phenomenological theory of superfluidity was constructed by L. Landau P[[in| in the 1940s. 
His key idea consisted in that He-II is a quantum fluid, in which motions can be only collective, 
i.e., in the form of quasiparticles. L. Landau formulated a criterion, according to which a 
quantum fluid is superfluid (SF), if its spectrum of quasiparticles (the curve of the energy 
e of quasiparticles as a function of momentum p) does not touch with the momentum axis, 
i.e., it has no points with e/p 0. L. Landau proposed to describe a gas of quasiparticles 
as the "normal component" of the fluid, and the "rest fluid" represents the SF-component. 
This yields the two-fluid hydrodynamics of He-II which allowed one to explain a number of 
properties of helium [TTl [T2] . But this theory is phenomenological and does not clarify the 
structure of quasiparticles. Therefore, the analysis within the two-fluid hydrodynamics can 
lead to the erroneous look at phenomena in some cases. 

A microscopic theory of He-II was started in a number of works. F. London assumed [13] 
that the superfluidity of helium can be related to the presence of a condensate in it. On 
the basis of this idea, L. Tisza p3] constructed a two- fluid hydrodynamics of He-II, in which 
the SF-component coincides with the condensate. The last assumption turned out false; a 
more correct two-fluid model was developed by L. Landau [U [10]. A. Bijl [15] proposed an 



approximate structure of the wave functions of the ground and first excited states of He-II: 
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A — the interaction parameter. The microtheory obtained the strongest stimulus due to the 
works by N. Bogohubov [17] and R. Feynman [181 HI]- In particular, the field-theoretic 
(FT) approach [16] allowed one to calculate the spectrum of quasiparticles (phonons) and 
the relative number of atoms of a condensate (below — "value of a condensate") for the 
system of interacting bosons under condition for the interaction strength to be small. The 
FT-approach involves the operator formalism and relevant methods. Then it was developed 
in a number of works, in particular, in [20] [23j. 

More significant results were obtained within quantum mechanical (QM) approaches al- 
lowing one to calculate, in various ways, the full A^-particle wave functions (WFs) of He-II 
such as those of the ground state (\I'o), the states with one quasiparticle = V'k^o, and 
the states with many quasiparticles 

n 

^(ri,...,r;v) = n(^kJ"'^o. (2) 

i=l 

It is worth noting that the interaction of quasiparticles in the structure was not considered 
([2]). The first sequential works in the QM-approach were executed by R. Feynman [T8l [19] . 
By using the variational method, he obtained the approximate formula 

for the spectrum of quasiparticles (SQ) of He-II, E{k), which indicated, for the first time, 
that the roton minimum of E{k) is related to a maximum of the structure factor S{k), i.e., 
to a short-range order in the fluid. 

The first works by R. Feynman were followed by the articles by N. Bogoliubov, D. Zubarev 
[T7] and R. Jastrow [2l] published simultaneously. In work [15], the logarithm of was 
sought in the form of the expansion of ([1]) in the interaction parameter. It was shown in [IT] 
that such a theory of perturbations is not fully correct. There, it was refined, and a new 
model of weakly nonideal Bose-gas combining the QM- and FT-approaches was proposed; 
the zero approximation ip-^^ = p j^ was deduced, and the SQ of a Bose-gas, 



M^^) +2Mfc)^, i^ik)= [V{r)e-^^dr, (4) 



coinciding with ([3]) at the weak interaction was found. Work [T7] induced the development of 

TV .1 

QM-approaches with the use of collective variables pj^ = J2 e~ J (r^ — coordinates of 

i=i 

helium atoms). However, that work did not else show that the logarithm of \E'o can be most 
explicitly expanded in the correlation series 

In^'o = Xl"2(ri - rfc) + J2 ^(ri - ffe, - r^, r^- - r^) + . . . (5) 

j<k j<k<l 



Similar, but yet not fully correct, structure of WF \I'o was proposed in [23], where was 
sought in reduced form M3(rj,rfc,r;) ~ f{rkj)f{rji)f{rki) and was not clearly indicated, that 



just the logarithm of \l/o should be expanded (it was shown in flbl flT]). R. Jastrow did not 
know, probably, work [15j, where structure was almost deduced. The proper structure 
([5]) was proposed independently and almost simultaneously by E. Feenberg [25], [26], C.-W. 
Woo [27J) T. Regge [28|, and I. Vakarchuk, I. Yukhnovskii [29]. Moreover, works [29] present 
a method of calculation of the fine structure of \E'o and ^'j^ [see formulas and ([7]) below]. 

From the historical viewpoint, it is interesting to note that work [16] was reported by N. 
Bogoliubov at the Landau's seminar, where it was criticized. One of the results of that work 
was SQ (jl]) which has a minimum in the region of finite momenta at a certain potential V{r). 
L. Landau noticed, without any doubts, this fact and submitted his famous work [TU] for 
publication already in a week, where he predicted correctly the spectrum of quasiparticles of 
He-II measured later on (Fig. 1). However, the work by Landau contained no mention about 
the Bogoliubov's model ... 




Fig. 1: Experimental spectrum of quasiparticles of He-II. 

I also note that the Soviet researchers in the USSR met difficulties as for the publication on 
the western journals published in English. The direct publications in the West were forbidden 
by the internal censorship, and they were possible mainly with delay through several journals 
published in Moscow and translated abroad. However, due to the hostility between the 
Bogoliubov's and Landau's schools, the representatives of the former had no way to the 
world scientific literature, since this way passed through Moscow journals. For this reason, 
the works by I. Yukhnovskii and I. Vakarchuk [2n]^plJ, in which one of the best modern 
QM-methods (see below) on the basis of work [T7j was developed, remained little known. 
One would expect that the period of totalitarian regimes terminates, all journals will be open 
for the authors from any country, and science will be aimed at only constructive purposes. 

3 Modern quantum-mechanical models of the microstructure of 
He-II 

The development of ideas advanced by London, Bijl, and mostly by Feynman and Bogoliubov, 
has led to the appearance of a lot of articles, from which four QM-methods should be selected: 
the methods of collective variables (CV) [28]-[32l[8], "correlated basis functions" (CBF) [33] 
[36], "hypernetted chains" (HNC) [37]-[l0], and "shadow wave functions" (SWF) [3l]-[ll]. 



The method of CV was developed in [30]. There, the WFs of the ground and first excited 
states, \l/o and are sought in the representation of the collective variables pk in the form 
of infinite series 

*o = ^^, So = Y^ — - — '5(ki + ... + ki)ai(ki,...,ki)pki---Pki, (6) 

V Vra i>2 '■■ ki,...,ki^O 

V'k = P-k+ 2!ViV ^ki-kP-ki+ ^j]^^ Pki+k2-kP-kiP-k2 + ---' (7) 

ki^o ki,k2^o 

where the l-th term describes the /-particle correlations; S — Kronecker delta, N — the total 
number of helium atoms. By performing the direct substitution of \E'o dHD and to the 

exact A^-particle Schrodinger equation, we can verify that \E'o and ipy^ satisfy it and can find 
an infinite chain of linked equations for the functions ai and bi. These equations are exact, 
but the whole chain cannot be solved analytically or numerically. Therefore, it is truncated, 
and only the first terms describing two-, three-, and sometimes four-particle correlations are 
retained. Such equations can be already solved numerically [H [321 USj. The method allows 
one to reproduce the experimental SQ of He-II with good accuracy, as well as the energy of 
the ground state (Eq) and the values of condensates, but with worse accuracy. The analysis 
can be carried out without fitting parameters (FPs), by using the structure factor S{k), or 
with FPs, by using the pairwise potential V{r). 

In the CBF method, the solution for ip-^^ is also sought in the form ([7]), and the functions 
bi and E{k) are determined by the Brillouin-Wigner theory of perturbations with regard 
for only, as usual, two first terms in ([7]). The method allows one also to calculate the SQ 
with good accuracy. The best results are obtained in [35] without FPs and in [36] with 
several FPs. In [33]-[35], the WF \E'o was not calculated. These works used many-particle 
distribution functions under certain assumptions and did not apply FPs. In [36], \E'o is sought 
in the coordinate representation ([5]) with regard for two first terms with the functions U2 and 
U3 which are chosen in a definite form with several FPs. In CBF models, the four-particle and 
higher corrections in the expansion of E{k) were not considered, though taking the subsequent 
corrections into account changes considerably the curve E{k) |35j . 

The HNC models are developed in [23 Eg] and [23 SO] • In [371 [35] , the WF ^0 is sought 
in the form similar to ([6]), and the functions ai are determined from the chain of equations 
obtained by means of variation of the total energy. The SQ is calculated by the Brillouin- 
Wigner theory of perturbations. The analysis can be performed without FPs: on the basis 
of the experimental curve S{k) or the potential (then S{k) is deduced). In [391 HO], ^0 
is determined in a similar way, but in representation ([5|) and with the use of a potential. 
There, the WF ipy^ and the SQ are calculated within a special method, which is equivalent 
to the calculation of ([7|). In HNC models, the quantities E{k), Eq, S{k), and the value 
of condensate can be calculated quite satisfactorily. As drawbacks, I indicate that the HNC 
method itself contains certain approximations, and only two first terms are considered in 
series ©-([7]). 

Of interest is the SWF method, whose idea consists approximately in the account of the 
delocalization of atoms by the introduction of special simple "shadow" factors in the structure 
of \E'o and \E'j^ in order to partially take higher many-particle corrections into account. The 
quantity \E'j^ is sought in the form [121 



^k = / • • -^^^ ■ ■ [-'^Y.l'^j ~ ~ Y^Mrji) +Us{sji))], (8) 

j j<i 

aq = J2 exp [iq.Sj + iqA{q) J2 i^j - ^dH^ji)], (9) 



where — shadow coordinates, Ur{r) = [b/r)"^ — the McMillan function, Us{r) = VAziz{di ■ 
r)/d2 — the renormalized Aziz potential, and C,m,b,di, and d2 — free parameters. The 
function X{sji) contains still two FPs; whereas A{q) and, in fact, Ug and Ur — free functions. 
The quantity \E'o is also described by formula ([8]) with aq = 1. The structure ([8]), IQ was 
not strictly deduced, but it is substantiated by a number of arguments with the use of the 
analogy with the formalism of path integrals. At certain values of FPs, it is possible to well 
reproduce the experimental data on E{k), Eq, S{k), and the value of a condensate. A part 
of FPs is determined under the modeling of \E'o) and the rest of FPs — in the calculation of 

\E'l^. Note that with 5"q = X^exp (iqSj), i.e., without use of FPs in the derivation of \E'jj., the 

j 

SQ is reproduced satisfactorily. In [43j, aq was replaced by (Tq_q'crq', which improves the 
agreement for E{k). It is probable that the result will be more exact, if aq will be represented 

by a series of the form ([7]) with o"q = X^^xp («qSj) instead of pj^. In our opinion, the SWF 

j 

method in such a form does not allow one to significantly improve the microtheory of He- 
ll, since structure ([S]), involves only partially higher correlations (as the first iteration 
^41] of the Schrodinger equation written in the form of a continual integral), whereas the 
complete consideration requires to deal with infinite series. In addition, the coefficients in ([8]) 
are not obtained, so that one has to use many free parameters and functions. This results 
in the obscurity whether the agreement with experiment is improved due to a more exact 
description of the microstructure of helium or simply due to a successful selection of a number 
of FPs. 

Thus, all models take only two- and three-particle corrections into account, by neglecting 
many-particle ones. In this aspect, the models are equivalent. The CV method seems to me 
the most transparent — it allows one to obtain the equations for the coefficient functions ai 
and hi, the required quantities E{k), Eq, S{k), and the values of condensates exactly and in 
a comparatively simple way. 

As an efficient numerical tool for the description of the microstructure of He-II, I mention 
the Monte-Carlo method (MC), in which the problem of next corrections does not arise, since 
the full solution of the Schrodinger equation is determined numerically at once, though only 
for the ground state |1H]^[S2]- The MC method gives the most exact description of the 
ground state, its energy Eq, the structure factor S{k), and the values of all condensates. But 
the method is not also perfect. In particular, it involves a potential V{r) which is not known 
exactly for small r. True, at a high barrier Vq = V{r 0), all potentials are close to the 
potential of a hard ball; therefore, the results of calculations for different (but high) barriers 
depend on the barrier height not so strongly. However, some arbitrariness concerning the 
potential can be noted: e.g., the Aziz potentials include several FPs, and one usually selects 
that Aziz potential among several ones which gives the better agreement of theory with 
experiment. Another drawback of the method consists in that it is numerical and does not 
allow one to see the analytic structure of a solution and its details. This yields the third 
drawback: due to the error of calculations, the method does not allow one to distinguish two 
cases where the value of a three-particle (for example) condensate is exactly zero or slightly 
different from zero, say, it equals 0.01. From the physical viewpoint, these cases are basically 
different. Moreover, the method does not yield the curve E{k) (and all characteristics of 
helium for T > 0), because the Schrodinger equation involves only the energy E, rather than 
the momentum k. Therefore, by minimizing E, one determines only E for k = 0. The whole 
curve E{k) must be sought involving the other methods, by choosing, for example, tp]^ in the 
form ([7]) with variational parameters. 

The hterature contains the discussions about the nature of a roton [IHl HH H US] . A roton 
is described by WF ([7]) which is the exact solution of the A^-particle Schrodinger equation; 
the same function ([7]) but with small k describes long-wave phonons. This implies that a 



roton is a phonon, but with short wavelength. To be more exact, a roton is a wave packet 
— a superposition of infinite in space waves ([7]) with various k in the interval Ak. This wave 

packet is compact, with a size of Ar ~ lOA [17], which yields a/c ~ 1/Ar ~ 0.1 A . That is, 
a roton is a collective excitation, rather than a one-particle excitation or a vortex ring. 

In many works, the potential V{r) of the pairwise interaction of He^ atoms is considered 
reliably determined. At large distances, r > 2.6 A [53], the potential is small < 10 K), 
negative, and well-known. But, at small r, the potential is positive, very large (Vq ~ 10^ -r- 
10^ i^'), and unknown exactly. Indeed, the electron shells of atoms overlap at small r, the 
interaction of electrons and the nuclei of both atoms is strong and is modeled [SI] in the self- 
consistent field approximation, which gives only an approximate estimate of the potential 
(with Vo ~ 10^ K), rather than the exact potential. The barrier Vq ~ 10^ K is also calculated 
[55j by approximating the potential V{r > 2 A) to the region of small r with the use of the 
quasiclassical approximation, whose exactness is insufficient for the given problem. In view 
of the approximate nature of the construction of the potential at small r, it is surprising that 
many works consider the potential at small r to be reliably established and to equal the Aziz 
potential ^53j. 

On the convergence of QM- methods. In series ©-([T]), the corrections ~ A^, and there 
is no small parameter in the equations for ai and bi. Nevertheless, the solutions for E{k), 
Eq, S{k) and the values of condensates should converge to the exact solution with regard for 
the increasingly greater number of correlative corrections in ©-([Tj). In the chain of linked 
equations for the functions ai [30j, each of the equations couples the functions ai-i, ai, and 
a;+i (and also nonlinear "lower" corrections). Despite the absence of a small parameter in the 
problem, the solution for the first functions (02, 03) should converge to the exact one if the 
increasingly greater number of equations in the chain is solved. This convergence is related 
not to a small parameter, but to the consideration of higher many-particle correlations. 

What does the theory allow one to calculate? In the frame of QM- methods, the SQ can 
be reliably calculated with regard for only two- and three-particle correlations in (jSD^dZ]) and 
without use of FPs [H [33]-[35l [38]. The deviation from experimental data ~ 10%. The 
quantity Eq and the value of a one-particle condensate ni are determined with significantly 
less accuracy, for which the experiment gives Eq = —7.16K, rii ^ 0.05 -7- 0.1. These values 
can be calculated better within the MC method. For the condensates, the problem consists 
in that values of the one- and two-particle condensates at T = are calculated [311 [321 [56] 
as the exponent of a sum containing all a;. Since the corresponding corrections are not small 
(~ 1/2), taking the following ai into account changes significantly the value of a condensate: 
on the transition from the two- to three-particle correlations, this value varies by several 
times [32l [56] . Thus, if one does not use FPs, then the value of a condensate in the QM- 
approaches (except the MC method) is calculated with low degree of accuracy. To the best of 
my knowledge, the FT-approaches do no contain also a quite correct scheme of calculations of 
the values of condensates for He-II. The value of a two-particle condensate n2 was calculated 
only in [32l [56l [57] , and no reliable experiment for the determination of n2 has been proposed 
till now. The structure factor S{k) is determined from the interaction potential, and all 
methods give a good agreement with experimental data. As for the potential, the situation 
looks more complicated. In the region of the potential well, r > 2.6 A, the models agree with 
one another, and the potential is satisfactorily reconstructed from the structure factor [S]. 
But the height of the potential barrier Vq is described by different models differently: from 
~ 100 ^ 1000 [211 [231 [581 [HI iS] to ~ lO'^K (with the Aziz potential [53]; in the frame 
of the MC and SWF methods, such a potential leads to a good agreement of theory and 
experiment for Eq, S{k) and the value of a condensate rii). The great difference of Vq, by 3-4 
orders, can be related to the efficiency of the potential at small r, the neglect of many-particle 
corrections in some models f2T[ [221 ESI [SI US] (which can cause the underestimation of the 



potential by 1-2 orders), or the fact that the Aziz potential is overestimated by 1-2 orders 
due to approximations made in its derivation. 

Let us summarize. The general problem concerning all QM-methods consists in that 
the exact solution is represented by the infinite series ©-([T]), for which all terms cannot be 
calculated. There are several ways to solve the problem: 1) To realize the break-down in the 
microtheory, by developing a basically new method of solution or the new form of a solution. 
For example, it is desirable "to convolve" series (l5])-(l7j) to simpler functions which would 
be exactly determined. But such a method is not discovered else. 2) To gradually calculate 
more and more corrections in (jSD^dZ]) without FPs or using them only for the potential to 
a minimal extent. The convergence of E{k), Eq, S{k) and the values of condensates can be 
clearly manifested already after the consideration of, say, 8-particle corrections. 3) However, 
a system of equations involving all corrections up to the 8-particle one is awkward and 
complicated. At present, in order to treat only the 5-particle correction, it is necessary to 
increase the operating rate and the random-access memory of modern personal computers by 

3 orders [15]. It seems that such great efforts are not expedient. Indeed, the studied two- and 
three-particle approximations give a satisfactory agreement between theory and experiment, 
and it is clear that an increase in the number of considered corrections will improve the 
agreement. 

Sometimes, the authors improve the agreement with experiment with the help of several 
FPs. Of course, using three and more fitting parameters, it is possible to "explain" any 
simple curve such as, e.g., the spectrum of quasiparticles of He-II. However, many-particle 
corrections in ©-([T]) are not small and, hence, are of importance for the solution. If they will 
be simply neglected, and the curve E{k) will "be constructed" with the help of several FPs, 
we will face with a successful fitting, rather than with a description of the real microstructure 
of He-II. In my opinion, to use more than one FP in a model is not expedient. 

The FT-methods [HI EDI ^[23], as far as I know, meet the same difficulty concerning He-II 
— the input exact equations following from the first principles cannot be solved exactly, and 
one is forced to omit many-particle corrections and to solve "truncated" equations. In this 
case, one fails to correctly introduce a small parameter in the problem (mainly due to the 
strong interaction of atoms), so that the neglected many-particle corrections are not small, 
generally speaking. 

I mention also the idea of He-II as a quantum crystal [59l [60] . The authors |60] consider 
that such disordered system as a fiuid cannot satisfy the requirement for a system to be 
completely ordered at T — (the third principle of thermodynamics). In my opinion, helium 
has no crystalline lattice. Otherwise, the SQ would be strongly anisotropic like that for the 
well-known crystalline phases of helium [^I]. In addition, the superfiuidity was observed in 
films of helium with thicknesses down to 0.1 f62] and 0.01 ^3] of the atomic layer at various 
mean distances between atoms, whereas this distance for a crystal must be of a single value 
equal to the lattice constant. 

4 Some unsolved problems 

At the present time, the structure of the composite condensate of helium is not quite clear in 
that how many atoms are in the one-, two-particle, and higher s-particle condensates (s > 3) 
and why the structure is just such. It is shown in work [32J in the three-particle approximation 
that the values of three-particle and higher condensates equal zero in the ground state of He- 
II. I note that the three-particle correlations were considered, so that it could be expected 
that the value of the three-particle condensate is nonzero. I presuppose that the values of 
three-particle and higher condensates turn out equal zero with regard for all corrections (i.e., 
for the exact solution) at any temperature. It is possible that the absence of the higher 



condensates is related, in a way, to the fact that the function ai (and bi) from (E]), (JTj) in each 
equation of the chain [30] is coupled only with two "adjacent" a^. This pair originates from 
the pairwise potential. 

Above, quasiparticles meant phonons (including rotons). However, in addition to phonons, 
helium should contain also microscopic vortex rings, as a heat excitations, with the radius 
r > sA, which was noted by L. Onsager [46j long ago. Such rings were registered in experi- 
ments with ions [471 EH EH]. The microscopic description of rings including the construction 
of the full iV-particle WF of a ring was not proposed till now. Only a solution in the mean- 
field approximation is available [66J. In this case, the A^-particle WF of helium with a single 

N 

ring is sought as the product H ^(r^j t) N of identical functions \E'(r, t) satisfying the Gross- 

Pitaevskii equation [67]. The exact parameters of the least ring (its radius, energy, and 
momentum) are unknown. Respectively, the contributions of rings to the heat capacity of 
He-II and to the A-transition are unknown as well. The difference of the experimental heat 
capacity of helium and that of the ensemble of phonons indicates that the contribution 
of rings to the heat capacity at temperatures T = 0.8 -h 1.2 K does not exceed 6 %. However, 
it can be significantly greater near T^. It is possible also that the rings appreciably affect the 
value of Tx- If this is true, one should observe an anomaly on the curve Tx as a function of 
the thickness c? of a thin film of helium at c? ~ (ic ~ 6 A [SH] {dc — the diameter of the least 
ring), which can be verified in experiments. 

The nature of the A-transition is not quite clear as well. A number of models were proposed, 
and some of them well describe the critical indices. But those models are different from the 
physical viewpoint. In particular, some authors believed that the rings play an important 
role in the A-transition (see references in [69]), but the other ones hold the opposite opinion. 
The question about the nature of the A-transition is fundamental, because it concerns, in 
fact, the nature of superfluidity. It is considered to be well established that the superfluidity 
is related to 1) all motions in He-II are collective, and 2) the spectrum of quasiparticles 
satisfies the Landau criterion. However, the superfluidity is already absent at a temperature 
somewhat greater that Tx- So that, we must answer the questions: What was changed in 
the medium on the microlevel, and why was the superfluidity broken? This is not quite 
clear. The process running at Tx is named sometimes the appearance of microturbulence. 
From the mathematical viewpoint, this means apparently that 1) structure ([2]) stops to be 
a solution of the A^-particle Schrodinger equation and 2) a long-range order (LRO) in the 
medium disappears. The LRO is determined for the s-particle density matrix Fg of He-II as 

lim F,(ri,...,r,|r;,...,r'J|B = F,(ri,...,r,)F;(r;,...,r'J, (10) 

\ri-r.\^co 

lim |F,(ri,...,r,)|2 = F,(oo) = [Fi(oo)]' = const > 

(in the 3D case), where i,j = 1, . . . ,s, and condition B means that |rj — rj| and |r'j — r^l 
are fixed. In a three-dimensional space, just the presence of the LRO causes the appearance 
of all condensates [311 [321 [56]. However, the things are not so simple with a LRO. Known 
are the theorems which forbid the existence of LRO in the 2D space in an infinite system. 
However, all systems in the nature are finite in all directions, and the values of a condensates 
and the form of LRO for such systems are unknown at present. 

This problem was considered for periodic boundary conditions (BCs) [TUl [7T] which are 
used explicitly or implicitly in all microscopic models of He-II. But since the periodic BCs 
are not realized, obviously, in practice for systems consisting of a great number of particles, 
the problem could be solved more exactly for zero BCs. In this case, the minimal momentum 
of a particle is not zero, but it is equal to ^tt/Lj {Lj — sizes of the system, L^, Ly, L^), 
which follows also from the uncertainty relations. The condensate is a fine effect, and we may 



expect that the proper setting of BCs is important for its evaluation. Apparently, both one- 
and two-particle condensates exist in finite systems, including finite 2D systems, and "are 
dispersed" over a number of the lowest levels. In the 3D space, the summation over all levels 
will result, apparently, in the condensate value -^j^^q for an infinite system. 

In 2D films, the superfluidity is conserved [62l[63]. If the condensates and LRO exist in such 
films, like the 3D case, we may conclude that both the condensate and long-range order are 
the microscopic peculiarities integral to the property of superfluidity. For bounded systems, 
the criterion for LRO flTUl) must contain |xj — Xj| L^^q < instead of |rj — r^l — > cxd, and 
the analogous substitution will hold for y and z. 

It was observed recently that He possesses electric properties [72]-[71], including those 
without external electromagnetic fields [72l [73] . For the time being, these interesting experi- 
ments have no explanation. The electric signals [72l[73] are likely related to the "spontaneous" 
polarization of helium which appears under certain conditions without external electromag- 
netic fields, as a result of the mutual polarization of helium atoms [751 176] . 

5 Conclusion 

I have tried to describe the state of the theory of superfiuid He-II. As seen, the theory is 
gradually developed, but many point remain obscure, including the basic properties of the 
system and the very property of superfluidity. This presents a wide field for the future studies 
and, of course, surprises. 

The work is performed under the financial support of the Division of Physics and 
Astronomy of the NAS of Ukraine in the frame of a Target program of fundamental studies 
N 0107U000396. 

^ Electronic address: mtomchenko@bitp.kiev.ua 

* A numerical expression of the term "several words" is not strictly definite in the literature, and I set it to be 
approximately 5000 words. 
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